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We revisit and extend previous calculations of glueball decay rates in the Sakai-Sugimoto 
model, a holographic top-down approach for QCD with chiral quarks based on D8-D8 probe 
branes in Witten’s holographic model of nonsupersymmetric Yang-Mills theory. The rates 
for decays into two pions, two vector mesons, four pions, and the strongly suppressed decay 
into four tt® are worked out quantitatively, using a range of the ’t Hooft coupling which 
closely reproduces the decay rate of p and w mesons and also leads to a gluon condensate 
consistent with QCD sum rule calculations. The lowest holographic glueball, which arises 
from a rather exotic polarization of gravitons in the supergravity background, turns out 
to have a significantly lower mass and larger width than the two widely discussed glueball 
candidates /o(1500) and /o(1710). The lowest nonexotic and predominantly dilatonic scalar 
mode, which has a mass of 1487 MeV in the Witten-Sakai-Sugimoto model, instead provides 
a narrow glueball state, and we conjecture that only this nonexotic mode should be identified 
with a scalar glueball component of /o(1500) or /o(1710). Moreover the decay pattern of 
the tensor glueball is determined, which is found to have a comparatively broad total width 
when its mass is adjusted to around or above 2 GeV. 
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I. INTRODUCTION 

The nonabelian nature of Quantum Chromodynamics (QCD)—the theory of the strong 
interactions—makes it possible to form bound states of gauge bosons, the so-called glueballs 
m- In pure Yang-Mills theory, these are in fact the only possible particle states and their spec¬ 
trum has been studied in detail in lattice gauge theory IIHS]. Glueballs are obtained for a range 
of quantum numbers , where J denotes total spin, P parity, and C charge conjugation; the 
lowest glueball has the quantum numbers of the vacuum, In the presence of quarks, 

the situation becomes complicated because glueballs can mix with qq states of the same quantum 
numbers. Lattice simulations of QCD including quarks are more difficult, but recent unquenched 
calculations continue to indicate the existence of glueballs [7] with the lightest glueball around 
1600-1800 MeV. 

The identification of glueballs in experimental data, however, remains elusive [HHII] and will 
be one of the objectives of the PANDA experiment at FAIR Experimentally, the O^"*" 

meson sector turns out to be particularly challenging. Listings of the Particle Data Group (PDG) 
[II] contain five isospin-zero scalar states in the energy region below 2 GeV: /o(500) or a, /o(980), 
/o(1370), /o(1500) and /o(1710), with the last two rather narrow states being frequently discussed 
as potential candidates for states with dominant glueball content |15H24j . Alternative scenarios 
with broad glueball resonances around 1 GeV and mixing with the a meson are also discussed in 
the literature [25l - [28] . 

A similarly unclear situation is found in the case of the lightest tensor glueball: lattice simula¬ 
tions obtain a mass between 2.3 GeV and 2.6 GeV |3HZ] while the PDG lists /2(1950), / 2 ( 2010 ), 
f 2 ( 2300 ) and /2(2340) as established states around and above 2 GeV, with several needing confir¬ 
mation [e.g., the narrow /j(2220) state that may have spin two or four, or may not exist at all]. 
Various approaches to low-energy QCD have also been applied to this region [29lf33] but a clear 
identification of a tensor glueball in the meson spectrum is missing. 

A central difficulty is the paucity of theoretical predictions of glueball couplings and decay 
rates from first principles. Lattice gauge theory provides information on Euclidean correlators and 
the extraction of real-time quantities is involved and fraught with uncertainties. Glueballs are 
particularly difficult to study when dynamical quarks are included. 

A completely different approach to strongly coupled gauge theories has been developed over 
the last one and a half decades in the form of anti-de Sitter/conformal field theory (AdS/CFT) 
correspondence, or more generally, gauge-string duality [MIES]. The AdS/GFT correspondence 
posits a map of correlation functions of gauge invariant composite operators with large number of 
colors Nc and large’t Hooft coupling to perturbations of certain backgrounds in classical (super-) 
gravity. Already in 1998, Witten [36] proposed a top-down construction of such a duality based on 
type-IIA supergravity, where both supersymmetry and conformal invariance are broken such that 
at low energies, below a Kaluza-Klein mass scale Mkk; the dual gauge theory is four-dimensional 
large-Vc Yang-Mills theory. The calculation of glueball spectra from type-IIA supergravity was in 
fact one of the first applications of “holographic QGD” [371 - 142] . (Glueballs have subsequently been 
studied further in more phenomenological, bottom-up holographic models in, e.g., Ref. [431 - 145] . 1 

Quarks in the fundamental representation can be added to the AdS/GFT correspondence in the 
form of probe flavor D-branes [46] . In type-IIA superstring theory there are D-branes of even spatial 
dimensionality, and the first attempt to include quarks in Witten’s model of nonsupersymmetric 
Yang-Mills theory was based on D6 branes m- This made it possible to study chiral symmetry 
breaking in the case of one flavor, which however did not permit a correct generalization to flavor 
number Nf > 1, an issue that was solved in 2004 by Sakai and Sugimoto [48l |49] by adding pairs of 
D8 and anti-D8 branes intersecting the color D4 branes of the Witten model. This model has been 
remarkably successful in reproducing various features of low-energy QGD while being firmly rooted 
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TABLE 1. Notations used here versus notation in Hashimoto et al. [SD], Brower et al. [55], and Sakai & 
Sugimoto [5S1155] 


in string theory with a minimal set of parameters—for given Nc and Nf, the only dimensionless 
parameter is the’t Hooft coupling A at the Kaluza-Klein scale Mkk- 

In this paper we shall use the Witten-Sakai-Sugimoto model to study glueball-meson interactions 
and to calculate glueball decay rates from the resulting effective interaction Lagrangians. This was 
first carried out by Hashimoto, Tan, and Terashima in Ref. [50], whose calculations we repeat (with 
important corrections) and extend. 

In addition to the lowest glueball mode in the Witten model, which happens to be rather 
different from the dilaton mode that plays this role in simpler bottom-up models of holographic 
QCD, we consider the (predominantly but not purely) dilatonic mode of the Witten model, as well 
as the tensor glueball and their excitations. We calculate decay rates into two and four pions, and 
we confirm the prediction of Ref. [50] that scalar glueball decay into four vr^ mesons is suppressed by 
evaluating the rate quantitatively. The latter receives contributions from multi-glueball interactions 
as well as from higher-order terms in the DBI action of the D8 branes, with the latter yielding the 
dominant piece. 

One of the main conclusion of our work is that the lowest gravitational mode in the Witten- 
Sakai-Sugimoto model appears to be ill suited to model the lowest glueball of QCD as found in 
lattice simulations, while the dilatonic mode has reasonable properties regarding its mass and 
decay rates. The lowest mode either has to be discarded on grounds of its exotic polarization along 
the compactified dimension of the type-IIA background or perhaps could find a physical role as a 
pure-glue component of the cj-meson [25] (which itself is absent in the Sakai-Sugimoto model) or 
the “red dragon” of Ref. [26] . 

We also make quantitative comparisons with experimental data on glueball candidates among 
scalar mesons at or above 1.5 GeV by extrapolating the mass of the holographic glueball and 
assuming weak mixing with qq states as the latter is parametrically suppressed at large W m 
and thus also in the Witten-Sakai-Sugimoto model m- Moreover, the decay pattern of the tensor 
glueball is worked out in detail, where also extrapolations to decays into massive pseudo-Goldstone 
bosons appear possible. 

In view of Refs. [241 152] . a particularly interesting feature of the holographic approach is that 
it admits narrow glueball states in the mass range predicted by lattice simulations, while the 
prediction of the gluon condensate is small, close to its standard SVZ value [53] . 


II. THE WITTEN MODEL OF NONSUPERSYMMETRIC YANG-MILLS THEORY 

The Witten model of nonsupersymmetric (and nonconformal) Yang-Mills theory in 3-|-l di¬ 
mensions is based on the AdS/CFT correspondence for a six-dimensional (0,2) superconformal 
field theory that is obtained from a large number Nc of coincident M5 branes in 11-dimensional 
M-theory. Their near-horizon 11-d supergravity geometry is the product space AdS 7 x S'^ with a 
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curvature radius L of the AdSy space that is twice the radius of the S^. With M5 branes extended 
along the directions 0,1,2,3,4, and 11, the line element of this space reads [S3] 

ds^ = ^ [r]f,t,dx^dx'' + {dx‘^)‘^ + + ^dr"^ + ^dQl, ( 2 . 1 ) 

where /i, = 0,..., 3 are (3+l)-dimensional indices (following [S3] we are skipping the index value 

10). The six-dimensional gauge theory living on the boundary of AdSy is a rather elusive maximally 
supersymmetric conformal field theory without a Lagrangian formulation. Dimensional reduction 
on a supersymmetry preserving circle with 

+ 27ri?ii, Rii=gsls, = a (2.2) 


leads to the near-horizon geometry of (nonconformal) D4 branes of type-IIA supergravity, whose 
dual theory is a five-dimensional super-Yang-Mills theory. 

Already in 1998, Witten proposed to use this correspondence as a basis for a holographic 
model of the low-energy regime of pure-glue Yang-Mills theory by a further circle compactihcation 
which breaks supersymmetry in the same way as supersymmetry is broken in the imaginary-time 
formulation of thermal field theory. The fermionic gluinos are subject to antiperiodic boundary 
conditions and thus become massive at tree level, whereas adjoint scalars acquire masses through 
loop corrections since they are not protected by gauge symmetry. In the limit of large Kaluza-Klein 
mass scale, the only remaining degrees of freedom are the gauge bosons. The dual geometry is 
given by a doubly Wick-rotated black hole in AdS 7 x 


ds^ = ^ [f{r)dxl + 'q^ydx^dx'' + dx\^ 


L2 

L2 dr2 L2 
j{r) 4 


'4) 


with /(r) = 1 - and a would-be thermal circle 


x^ ~ x^ -|- 27ri?4, Ri = 


1 


L2 


Mkk 3rKK ’ 


(2.3) 


(2.4) 


where the relation between xkk and Mkk is determined by the absence of a conical singularity at 
'T = xkk- The background also has a Ramond-Ramond (R-R) nonvanishing antisymmetric tensor 
gauge field with units of flux through the S^. 

The relation to the type IIA string-frame metric is 


ds^ = Gj^jydx^dx^ = e guNdx^dx^ -|- e^^^^{dx^^ + Aj^dx^)'^, (2.5) 


with M, N = 0,... 9 and M, N additionally including the index 11. This leads to a nonconstant 
dilaton = (r and Am = 0 for the above background geometry. 

For later use we introduce the alternative radial coordinates U G {Ukk,oo) and Z G (0,oo), 
used also in Refs. [ISIIIH])) through 


U = 


r2/2 

L ’ 


K{Z) = 1 + Z^ = 


KK 


TT3 

'^KK 


( 2 . 6 ) 


Note that the holographic boundary is at infinite values of r, U, and Z. 
In terms of the radial coordinate U the 10-dimensional metric reads 


ds^ = 


U \ 


3/2 


—J [7?^^dx^dx^ + /([/)(dx^)2] + y-jpj + U^dnl 


RD4: \ 


3/2 


dU^ 


(2.7) 
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with f{U) = 1 — (Ukk/U)^', the nonconstant dilaton is given by 

( 2 . 8 ) 

The parameters of the dual field theory are given by [UHISI |55] 0 

9ym — 2t^R ~ (.^/2)^ = .Rd4 = '^9s]^Js- (2-9) 

At scales much larger than Mkkj the dual theory turns into 5-dimensional super-Yang-Mills the¬ 
ory. However, it is not possible to make Mkk arbitrarily large without leaving the supergravity 
approximation. 

The dual gauge theory exhibits confinement. Wilson loops connecting heavy quarks at the 
boundary with large spatial separation along x are represented by fundamental strings that mini¬ 
mize their energy by having most of their length at minimal radial coordinate. The effective string 
tension therefore tends to the value 


a = 



2^12 9tt9xx 


U=Ukk 


1 (UkkV^^ 
2nlj { R J 


25ym-^c „ ,2 


( 2 . 10 ) 


In accordance with confinement, the dual theory has a mass gap for fluctuations of the background 
geometry with scale set by Mkk- 


A. Holographic glueball spectrum 

Ignoring all Kaluza-Klein modes on the compactification circles and all nontrivial harmonics 
on the with nonzero R charge, the bosonic normal modes of the supergravity multiplet can be 
interpreted as glueballs in the dual 3-|-1-dimensional Yang-Mills theory |371l42j 0 There are in total 
six independent wave equations for various scalar, vector, and tensor modes, which were denoted 
as S 4 , T 4 , V 4 , N 4 , M 4 , and L 4 in [l 2 ], see Table [TT} These give three distinct possibilities to obtain 
modes with = O"*"'*' quantum numbers, corresponding to the 3-|-1-dimensional scalars Gnqi, 
G 4 , 4 , and the 5^ volume fluctuation G“q,, where the index a refers to the S'^. The latter, termed 
L 4 in Table [III has a lowest mass eigenvalue ~ 3.57Mkk which is larger than those of all the other 
wave equations and will be ignored in what follows. 

The remaining two towers of scalar modes are described by the wave equations denoted S 4 
and T 4 . The lowest mass eigenvalue is found in S 4 , which corresponds asymptotically to 11- 
dimensionally traceless metric fluctuations in Gu, Cnqi, and G 44 . The other scalar mode does 
not involve G 44 and can be attributed to the dilaton derived from Gnqi. It is degenerate with the 
2~^~^ tensor mode (wave equation T 4 ) that is provided by transverse-traceless fluctuations in Gij, 
i,j = 1, 2, 3. (It is also degenerate with the vector mode l^^ derived from Gn^i, but this mode is 
discarded as spurious from the point of view of the 3-|-1-dimensional Yang-Mills theory because of 
negative “r-parity” [l 2 ], implying that its dual operator is odd under a reflection —)• —x^.) 

Pseudoscalar (0 '") modes are obtained from the 1-form field component 6*4 descending from 
Gii ,4 (wave equation V 4 ), whereas the 3-form field of 11-dimensional supergravity is responsible 
for vector modes: a vector 1^ from the antisymmetric tensor field Bij (wave equation N 4 ), and 

^ This is based on a normalization of the Yang-Mills action as — Mwhich differs, however, from the 
convention used in particle physics, where the coupling constant of SU(Ac) gauge theories is invariably defined 
as H so that = 2 (;ym- This means that the QCD coupling is given by as = g'^ji/l'K) — 

(/YM/(27r) = \I{2'kNs) in terms of the’t Hooft coupling A = N^g^ as used here. Since we do not attempt to match 
with perturbative QCD here, this is of no concern for the calculations performed below (it is, however, important 
to take into account when comparing quantitatively with weak-coupling results, see also footnote 1 in Ref. [56]). 

^ In Ref. [^ this analysis was recently extended to modes obtained by breaking the symmetry of the 5"^. 
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FIG. 1. The glueball spectrum of the Witten model (a) in units of Mkk, (“exotic” scalar modes in green), 
compared to the spectrum obtained in the recent large-fVc lattice calculations of Ref. [SH] (b) in units of the 
square root of the string tension yTr, juxtaposed such that the lowest tensor mode is matched. The dotted 
lines in figure (b) give the glueball spectrum of the Witten model when expressed in terms of the string 
tension of the Witten model with the standard set of parameters (3.8) for the Sakai-Sugimoto model. 


a vector 1 from the 3-form field components Cjj 4 (wave equation M 4 ). All other modes can be 
discarded due to negative r-parity. 

The glueball mass spectrum resulting from the numerical results listed in Table |TI] is displayed 
in Fig. where it is compared with recent lattice results at large Nc from Ref. [58| . which is in 
fact rather similar to that obtained for Nc = 3 mis]. When juxtaposed such that the lowest tensor 
mode is matched, the holographic spectrum roughly reproduces the pattern obtained in lattice 
gauge theory. Missing states of spin 2 with PC ^ and higher spin states might be due to 
closed string modes. On the other hand, there is a certain proliferation of O"’"*' states due to the 
existence of modes involving G 44 , which have been termed “exotic” in Ref. im, where they were 
first considered. In fact. Ref. mi suspected that only one of the towers of scalar states may survive 
in the limit Mkk 00 , where the Witten model would turn into an exact string-gauge dual of 
large-W Yang-Mills theory. 


B. Normalization of glueball modes 

In order to be able to derive effective actions of the glueball modes and their interactions, we 
need to calculate the normalization factors required for a canonical kinetic term. For this purpose 
it is convenient to use the 11 -dimensional notation, where the fluctuations take their simplest form. 
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Mode 

S4 

0 ++ 

T4 

0 ++/ 2 ++ 

V 4 N 4 M 4 

0-+ 1+- 1— 

L 4 

0 ++ 

n =0 

7.30835 

22.0966 

31.9853 53.3758 83.0449 

115.002 

n=l 

46.9855 

55.5833 

72.4793 109.446 143.581 

189.632 

n =2 

94.4816 

102.452 

126.144 177.231 217.397 277.283 

n=3 

154.963 

162.699 

193.133 257.959 304.531 

378.099 

n=4 

228.709 

236.328 

273.482 351.895 405.011 492.171 


TABLE 11. Our results for the mass spectrum of AdSy black hole metric fluctuations in the notation of 
[42] (i.e. in units of = -^kk/ 9) obtained by spectral methods cross-checked with a shooting method. 

The results for the lowest modes agree completely with Ref. [H], while for certain higher modes there are 
deviations in the last few digits. assignments are given only for the modes with even “r-parity” that 

are expected to have a counterpart in QCD. 


1. Lowest (exotic) scalar glueball 

The lowest scalar glueball is associated with fluctuations involving asymptotically (for 
r —>■ oo) (iG 44 = —45Gii = —4(5^22 = —4(5^33 = —45Giipi. In the bulk, other metric components 
are also involved, leading to the following “exotic polarization” m 


5Gaa = -jj^fHE{r)GE{x) 

2 

^G^i, = Y^HE{r 

5Gu,ii = '^^HEir)GE{x), 




0^6 


4 ' 5r6-2r|j^ 


f) f) 

M 


sa„ = -T/ 


-1 


34k 


5^6 _ 

sr = 90 

M|L2(5r6 - 2 r|j^j 2 

where the eigenvalue equation is given by 




HE{r)GEix), 


HEir)d^GEix), 


432 r2 


TL^Ml ) HE{r) = 0. 


,(5r6-2r6j^)2 

Integration over the 5^ reduces the 11-dimensional supergravity action to 

5 = ^{L/2)^nA j (fxyj- det G {r{G) + 

with 2^2^ = (27r)®/®5^ and 114 = 87r2/3. 

Inserting the metric fluctnations (2.11) into the 7-dimensional action gives 


d’xyj- det G {r{G) + 




with 


= -Ce / d^x dx^ - [{d^GE)'^ + M|G|] 


f^drr^B ^ 
Ce = I • 


( 2 . 11 ) 


( 2 . 12 ) 


(2.13) 


(2.14) 


rKK 


(2.15) 
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For the lowest eigenmode He we obtain numerically 

4 

Ce = 0.057395 [//^(rKK)]^^. (2.16) 

[This deviates from the result given in Ref. [50] by a factor ^ that seems to be missing in their 
Eq. (2.19).] 

Requiring that upon integration over and the scalar field Ge{x) is canonically normalized 
leads to 

[HE{rKK)]-^ = [He{Z = Q)]-^ = ^0.0097839 at, Mkk 

v 2 

= 0.0069183 A^/^iVcMKK- (2.17) 

(This differs from |50] only by the explicitly written factor 11^/2.) 


2. Scalar and tensor modes from the tensor multiplet 
A scalar mode O^"*" that does not involve metric components with index 4 is obtained fronij^ 


<5Giiqi = -2>ji^HE{r)G d{x) 


6 G 






dpdu' ^ . 

Gd{x) 


(2.18) 


Since upon reduction to 10 dimensions is essentially the dilaton, we shall refer to this mode 

as predominantly dilatonic. [Note that also the “exotic” mode (2.11) involves a dilaton component, 
but that there the dominant component is 5 G 44 . It should also be kept in mind that the attribute 
“exotic” only refers to the holographic origin of this mode, and not to any exotic quantum 
numbers in the dual field theory.] 

The tensor glueball 2++ is dual to metric fluctuations that have neither 5 G 44 nor (5Giiqi, but 
contain a transverse traceless polarization tensor in 5G^y. For example, one can choose as only 
nonvanishing components 


= — 6 G 22 = — j 2 ^T{r)GT{x). 

The radial functions Hd,t are determined by the equation 

- 4K)|:^D,r(r) + L^M^HD,T{r) = 0, 

with = M^) = My. 


Calculating the normalization of these glueball modes in analogy to (2.14) leads to 

dr ( 6HE{r)‘^ 


poo 

Cd,t = / 

JrKK 


L3 \ Ht{i 


For the lowest eigenmode He we obtain numerically 

Ct = 0.22547 [i7T(rKK)]'^ 

and an analogous result for Ce with a coefficient 6 times as large. 
This leads to 


[HD,TirKK)] ^ = [HD,TiZ = 0)] ^ = A^/^ A^cMkk I 


0.033588 

0.013712 


(2.19) 


( 2 . 20 ) 


( 2 . 21 ) 


( 2 . 22 ) 


(2.23) 


® As discussed recently in Ref. [SZl, more possibilities for scalar (and other) glueball modes are obtained if Ramond- 
Ramond field fluctuations which partially break the SO(5) symmetry are included. 
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C. Glueball field/operator correspondence 

The above metric perturbations are sourced by operators in the dual field theory, which is 
five-dimensional super-Yang-Mills theory compactified on the circle along 

The operator dual to the tensor perturbations is simply the five-dimensional energy-momentum 
tensor with three-dimensional indices. Omitting the adjoint scalars of the five-dimensional theory, 
we have 

= T™ + - \5mnFi^,F^^^ + ..., (2.24) 

where A 4 is a further scalar that like the adjoint scalars of the five-dimensional theory becomes 
massive through loop corrections. 

The operators dual to the exotic and the predominantly dilatonic scalar modes can be inferred 
from their couplings to the fields in the DBI action of D4 branes in the limit of r —?• 00 [^. The 
exotic scalar mode and the dilatonic one, turn out to source, respectivelylu 

+ F^^F^ - \fI + ..., (2.25) 

= +\f,,F^^^ - Ito™ + - \fI + .... (2.26) 

The difference = F^yF^^'^ is the purely four-dimensional glueball operator, which is 

dual to \5G^^ — 6G^jq. However this linear combination is not a normal mode in the gravitational 
background. We therefore need to keep the exotic and the predominantly dilatonic mode, of which 
both, or perhaps only one of them, might correspond to the glueballs of the four-dimensional 
Yang-Mills theory. To really end up with the latter, one would however need to take the limit of 
large Kaluza-Klein mass Mkk, which is necessarily leaving the supergravity approximation. In this 
limit, both modes will presumably receive important corrections. If one of the modes drops out of 
the spectrum, one might suspect that it will more likely be SGfjp^ as it includes a then spurious 
polarization component < 5 ^ 44 . 

In the following we shall consider both modes, as well as the tensor mode, when calculating 
glueball-meson interactions within the Witten-Sakai-Sugimoto model, extending the analysis of 
Ref. [50], which only studied the lowest (exotic) mode. 


III. THE WITTEN-SAKAI-SUGIMOTO MODEL 

Sakai and Sugimoto introduced chiral quarks in Witten’s model of pure-glue Yang-Mills theory 
by means of Nf probe D8 and anti-D8 branes that fill all spatial directions except the Kaluza- 
Klein circle |48[ I49| . Quarks and antiquarks are thus localized on separate points of the 4-|-l- 
dimensional boundary theory. The global flavor symmetry \]{Nf)L x \J{Nf)ii is however broken 
spontaneously, because the subspace x'^-U has the topology of a cigar forcing the D8 and anti-D8 
branes to join in the bulk. The action of the joined D8 branes which describes the dynamics of qq 
mesons through flavor gauge fields on the branes reads 

S'ds = — TbsTr J (f xe~^ — det {qmn + ‘^’^c^'Fmn) + Scs 

= -(2W)2rD8Tr I (fxe-'^^g (^1 + FpqFrs + 0(F^)^ + Kcs (3.1) 


^ We disagree here with Ref. |42| which attributed to SG^ and Too to SG^. 
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with Tds = (2vr)“®/7®, guN the metric on the 8+1-dimensional world volume induced by (2.7), 
and shifted such that e'*’ = Because no backreaction of the D8 branes on the 

10-dimensional background of the Witten model is taken into account, this corresponds to the 
quenched approximation of QCD, as indeed appropriate for the large-A^c limit at fixed Nj. (For 
attempts to go beyond the quenched approximation see Refs. [5^IHH] .! 

In the original version of the Sakai-Sugimoto model that we shall use here, the D8 and anti-D8 
branes are put at antipodal points so that they join at the minimal value U = Ukk- In this case 
it is most convenient to use the dimensionless coordinate Z = introduced already 

above, but extended to the range —oo ... + 00 so that the radial integrations of the D8 and the 
anti-D8 branes are combined. The part of the DBI action quadratic in the flavor field strength 
then reads 


SM’ = -'^Tr 




dZ 


-K-^l^rfPrT^F^^Fp^ 


+ KMl^rf'^F^zF.z 


with K = 1 + Z"^ and 




where (2.9) as well as = 87r^/3 and = (3/2)^17kk/Rd 4 have been used. 
The Goldstone bosons of chiral symmetry breaking appear as 


Sbs = 4- [ d^xTi {U-^d^uf + ..., U = Fexp\i 


dZAz 


which determines the so-called pion decay constant in terms of A and Mkk as 


f = 

J TT 




AW-TIkk- 


(3.2) 


(3.3) 


(3.4) 


(3.5) 


Massive vector and axial vector mesons arise as even and odd eigenmodes of A^if^ = i/jn{Z)vl^^ (x) 
with eigenvalue equation 

- (1 + ((1 + Z^)9zV'n) = AnV’n, V’n(±oo) = 0. (3.6) 

The lowest mode is interpreted as the isotriplet p meson (or the uj meson for the U(l) generator) 
with mass m?p = AiMkk with the numerical result Ai = 0.669314.... 

( 2 ) 

The next-highest mode vli with eigenvalue A 2 ~ 1.569 is an axial vector that can be identified 
|48| with the meson ai(1260). The experimental value for the ratio mai/mp ~ 1.59 is remarkably 
close to v^"^ 2 /Ai) ~ 1.53. Also the experimental value for the mass of the excited /9(1450) with 
mp*/mp K, 1.89 is is close to y^As/Ai) ss 2.07. This nice agreement may however be a bit 
fortuitous, since recent lattice simulations m at large Nc, extrapolated to zero quark mass, give 
the higher values ma^/nip k, 1.86 and mp*/mp ^ 2.40. This would correspond to errors 21% and 
16%, respectively, which may still be considered a success given that already the mass of v)i is 
above Mkk- (For more checks of the quantitative predictions of the Witten-Sakai-Sugimoto model 
see Ref. [62].) Optimistically, one can therefore hope that the Witten-Sakai-Sugimoto model is a 
useful approximation to QCD up to masses of two or three times Mkk- 
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A. Choice of parameters 


Matching the result for the p meson mass with its experimental value, rup = \/^Mkk ~ 
776 MeV Infixes the Kaluza-Klein mass to [38lll9] Mkk = 949 MeV. This determines the masses 
of the other vector and axial vector mesons, which come out in rough agreement with experiment. 
The masses of the lowest (exotic) and the predominantly dilatonic scalar glueball, the tensor 
glueball (degenerate with the dilatonic scalar), and the lowest pseudoscalar glueball are fixed to, 
respectively. 

Me = V7.30834/9Mkk ~ 855 MeV, 

Md = Mt = \/22.0966/9Mkk « 1487 MeV, 

Mp = ^31.9853/9 Mkk ~ 1789 MeV, 

Me* = \/46.9855/9Mkk « 2168 MeV, 

Me* = Mt* = ^55.5833/9 Mkk ~ 2358 MeV, (3.7) 


where we have also given the masses of some of the corresponding excited states (marked by a 
star). 


The lowest scalar glueball involving the exotic polarization (2.11) with a dominant < 5^44 com¬ 


ponent is found to be only 10% heavier than the p meson. This is in stark contrast to lattice results 
both for quenched Nc = 3 and Nc = 00 QCD [58], where the lightest glueball is about twice as 
heavy. 

A possible modification of the Sakai-Sugimoto model consists of choosing a nonmaximal sepa¬ 
ration of the D8-D8 branes [631164] . The latter then join at a value U = Uq > Ukk and the mass of 
a string stretched between Ukk and Uq has been interpreted as a “constituent” quark mass. Un¬ 
fortunately, this only makes the problem worse: Nonmaximal separation increases the eigenvalue 
Ai |65| while the glueball spectrum is unaffected. With a constituent quark mass of 310 MeV and 
keeping the mass of the p meson fixed as done in Ref. [66 |, Mkk is reduced to 720 MeV, which 


reduces all values in (3.7) by 25%. 


With maximal separation and the standard choice Mkk = 949 MeV, the mass of the dilatonic 
glueball is not far from the numerical result obtained in lattice gauge theory for the lightest scalar 
glueball state, while a degeneracy with the tensor glueball is not observed there—the latter is 
instead significantly heavier. This degeneracy might perhaps be lifted by higher-derivative cor¬ 
rections when going beyond the leading supergravity approximation. Similarly, it is conceivable 
that only the dilatonic glueball survives in the (unfortunately inaccessible) limit to a complete 
holographic QCD and that therefore the lowest scalar mode is to be discarded. We shall come 
back to this question when calculating the decay width of the various glueball states. 

In order to calculate glueball-meson interactions, we shall need to extrapolate to finite coupling 
and finite W = 3. The original [481 149] and most widely used choice is obtained from matching 
/jT ~ 92.4MeV in (3.5) which gives 


K = AiVc/(2167r^) = 7.45 • 10 


-3 


A PS 16.63 (Ac = 3). 


(3.8) 


[The original and published version of Ref. [48114^ contained an error in the prefactor of the D8 
brane action for Nf > 1 involving a different definition of k, which led to a’t Hooft coupling of 
about 8.3 and effectively a correspondingly reduced pion decay constant. This error, which was 
later corrected in the e-print versions of Ref. [481 149] , did not affect the mass spectra of mesons 
obtained in Ref. |48l|49|, but it does affect all interactions. Unfortunately, Ref. |50| still employed 


The mass of the lo meson, which is degenerate with the p meson in the Sakai-Sugimoto model, is only slightly 
higher in real QCD. 


5 
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the incorrectly matched’t Hooft coupling, affecting all meson and glueball decay rates calculated 
therein.] 

In what follows, we shall take (3.8) as the standard choice, but also consider as an alternative a 
value of the’t Hooft coupling obtained by matching mp/y^, where a is the string tension (2.10), to 
the large-A^c lattice result of Ref. [6lj. Ref. [61] obtained nip/^ = 1.504(50), whose central value 
corresponds to A = 12.55. With the “standard” value A « 16.63 the Sakai-Sugimoto model predicts 
vriply/a K, 1.306, which agrees within 15% but points to a smaller ’t Hooft coupling and thus a 
smaller string tension. A smaller’t Hooft coupling has also been argued for in Ref. ISZI, where the 
spectrum of higher-spin mesons obtained from massive open string modes has been considered. We 
shall therefore consider a downward variation of A ~ 16.63... 12.55 to get an idea of the variability 
of the predictions of the Witten-Sakai-Sugimoto model. 

Before turning to decay rates, we consider two other predictions of the Witten-Sakai-Sugimoto 
model at finite W where the concrete value of A matters. 

At inhnite Nc, the Goldstone bosons include also a massless rj' pseudoscalar meson from the 
spontaneous breaking of the Ua( 1) symmetry, whose anomaly is suppressed at W —)> oo. However, 
at finite the Sakai-Sugimoto model predicts a finite mass for the rj' meson through a Witten- 
Veneziano formula evaluated already in |48j with the result 


AMkk- (3.9) 

With Mkk = 949 MeV and A 16.63 (or 12.55) the numerical value for Nc = Nj = 3 turns out to 
be 967 MeV (730 MeV). The higher value is surprisingly close to the experimental value 958 MeV, 
but actually a smaller value than that might perhaps be expected given the absence of a strange 
quark mass. At any rate, the right ballpark seems to be reached with the parameters considered 
here. 

Another quantity of interest, in particular in connection with glueball physics, is the gluon 
condensate which was calculated in Ref. |55| as 

(3.10) 

For A Ri 16.63 this yields = 0.0126 GeV^, almost identical to the standard SVZ sum rule value 
|53j . while for A = 12.55 a significantly smaller value of 0.0072 GeV^ is obtained. Using sum rules 
both smaller [68| and larger [^2] values than the standard SVZ sum rule value are discussed in the 
literature, while lattice simulations typically give significantly larger values, which are however of 
the same size as ambiguities from the subtraction procedure m- While a quantitative comparison 
thus does not seem to be in order, we note that the gluon condensate is predicted to be small. 


rriri' = 



B. Normalization of qq modes 

For the calculation of decay rates we will initially consider Nf = 2, dropping the strange quark 
whose nonnegligible mass cannot be easily accommodated within the Sakai-Sugimoto model (see 
however Ref. mHH]) ; the possible effects of the finite quark masses will be discussed in Section [V} 
In the chiral Sakai-Sugimoto model, the Goldstone bosons are the massless pions contained in 

Tlz = UKKMZ)7r{xn, (3.11) 

where Ukk has been included to render the mode function 4>o{Z) dimensionless]^ The U(l) part 

® For our purposes it is most convenient to keep Az nonzero. The frequently adopted gauge choice Az = 0 leads to 
a different but physically equivalent field parametrization of the Goldstone bosons. 
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of Az corresponds to the r]' meson, which is a Goldstone boson only at infinite Nc\ for finite Nc it 
receives a mass through the Witten-Veneziano mechanism [48j (see Eq. (3.9) below). 

The only vector mesons that we shall consider will be the isotriplet p meson described by the 
traceless part of 




(3.12) 


and the isosinglet a; meson given by the corresponding expression proportional to the unit matrix. 

Following Ref. [50] (which here differs from |48(l49j) we choose the generators of the SU(2) flavor 
group such that TrT“T^ = <5“^. Canonical normalization of the fields 7r“ and in (3.2) such that 
upon integration over Z one has 


5 = -Tr / (Tx 


+ 2^1^KKpI, + • ■ 


leads to 


2k dZK-^l^(pl)xf = \, 


2K{UKKAlKKy 


dZK{<Pof = 1. 


(3.13) 

(3.14) 

(3.15) 


The first relation determines the value of at Z = 0 with the help of the numerical result 

poo 

/ dZ = 2.80302... (3.16) 


while the second fixes the normalization of (/>o oc 1/K = 1/(1 -|- Z^) as 

UkkMkk (po = 


\/ 2 'kk K 


(3.17) 


C. p and u) meson decay 


The p-TT interactions are determined by the second term of (3.2), usin^T)j^ = d^Az — dzA^ — 
i[A^,Az]. The effective vertex between the four-dimensional fields p and tt’s are obtained upon 
integration of the resulting products of the mode functions ipi{Z) and (po oc 1/K. For the process 
p —7- TTTr we need specifically 


C-p-K-K = -gp-K-K(-abc{dp,-n°‘)p^^T:^, Qpnn = ^ J = ^2 X 24.030 A 2Nc (3.18) 

This agrees with the numerical value given in table 3-34 of Ref. [l9| for = gp-ww- {gp-n-ir/V^ 

was denoted as cg in |50j : we will reserve Cj, f = l,2,3..., for the coefficients in the interactions of 
the glueball field with mesons, for which we will follow the conventions chosen in m-) 

The amplitude for the decay of a p meson at rest with polarization = (0, e) into two pions 
with momenta = (|p|,p) and = (|p|, — p) reads 


Ad. — igp-KTT Qp) - ‘^igp-KTT ® ■ P- 


(3.19) 


^ Here we follow the conventions of Ref. [50] . Note that in Ref. [481149) the matrix-valued flavor gauge fields are 
antihermitean. 
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The expression for the decay rate involves a directional average, leading to 

r / f -70 l-^l^ _ ^ / 0-1535 (A = 16.63) 

dvr J 167rm2 dSvr ~ XN^ ~ \ 0.2034 (A = 12.55) 


(3.20) 


which compares remarkably well with the current experimental value Tpjmp = 0.191(1) from 
Ref. [13] (although it should be noted that in this process the finite pion mass implies a reduction 
by about 20% compared to a decay into massless particles so that the coupling appears 
somewhat underestimated with our range of parameters for the Sakai-Sugimoto model). 

The decay of the ui meson into and which is due to the Chern-Simons part of the 

D8 brane action, has been calculated in |39|, with the result 2.58 MeV for the dominant 3-pion 
decay, which is significantly below the experimental value ~ 7.6 MeV. However, the result of |49j is 
proportional to A“^. Varying again A from 16.63 to 12.55 gives the range 2.58... 7.96 MeV, which 
happens to include the experimental value. 

So the model appears to make reasonable semi-quantitative estimates for meson interactions, 
which is quite remarkable given that after fixing the mass scale and setting Nc = 3, there is only 
one free parameter, namely A. This certainly makes it interesting to consider the predictions of 
this model for glueball decay rates in detail. 


IV. GLUEBALL-MESON INTERACTIONS 


The glueball modes, which have been obtained in Sect. IIA in terms of 11-dimensional metric 
perturbations translate to perturbations of the type-IIA string metric qmn and the dilaton 

$ according to (2.5). Explicitly, this gives 

-..3 r 


9fj,u — 


L3 


L2 


L2 


1 + ^'^‘^1141 j VtMU + 


r^f f L? L? 

544 - 

L / r'^ f 

9rr = ^ (^1 + 


9r^ — 


9fin = 


L 




m$/3 _ 

L2 


1 + 


(4.1) 


Here we differ from Ref. |50| where the metric fluctuations g^iQ on the S‘^ have been omitted. As 
one can check (Appendix]^, the 10-dimensional equations for the glueball modes are satisfied only 


when the fluctuation in is keptjZ] 


We shall consider in turn the lowest glueball dual to the metric fluctuations (2.11), referred 
to as “exotic” because it involves 6 G 44 besides dilaton fluctuations in JGnqi, the predominantly 


dilatonic glueball associated to (2.18), and the tensor glueball with metric fluctuations (2.19). 


Inserting the respective metric fluctuations in the D8 brane action and integrating over the bulk 
coordinates yields effective interaction Lagrangians which are given in full detail in Appendix [B| 


In 10 dimensions, the induced fluctuations in gnn are in fact necessary to decouple the mode L 4 , which in 11 
dimensions corresponds to pure 5"* volume fluctuations, as can be seen from the explicit 10-dimensional calculations 
in Ref. [5^ . 
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p a 


q b 


FIG. 2. Leading-order glueball decay into two pions. 



vertex 

value 

CIIV2 

GEdirdir 

44.304 A -5 

C 2 I V2 

GePP 

5.0318 A -5 

C 3 / 

GEdpdp 

49.334 A -5 A-iMj^k 

C4/-\/2 

GEpdp 

-7.4810 A -5 A-iM+j^ 

C 5 / \/2 

GEPT^dlT 

1428.1 A-i AT^Mj^k 

ci/y2 

GEdirdir 

11.590 A-1 A-iMj^k 


GePP 

2.0970 A-5 A-iMj^k 

c-ij\f2 

GEdpdp 

12.814 A -5 A-iMj^k 

C 5 / 

GEPT^d-K 

359.33 A-i AT^Mj^k 

cI/V2 

G^dirdi: 

24.641 A -5 A-^Mj^k 

cl|^f2 

G*ePP 

-0.8227 A-5 A-iMj^k 

C 3 /V 2 

G*Edpdp 

27.906 A-5 A-iMj^k 

c\l^ 

G*Epdp 

-1.7467 A-1 A-iM+j^ 

cH^ 

G*EPTrdTr 

858.66 A-i AT®M-^ 


G*EdTTdTT 

4.5843 A -5 A-im-^ 


G*ePP 

-1.2390 A -5 A-iM-^ 

CI/V2 

G*Edpdp 

5.3829 A -5 A-^Mj^k 

cl/^ 

G*EPEdTT 

176.99 A-i AT^Af-^ 


TABLE III. Coupling coefficients in interaction Lagrangian of lowest glueball. (Here we give numerical 
results for c ,;/\/ 2 to permit a comparison with the results fisted in m, with which we disagree by a factor 
of in (2.17). Taking this into account we agree with all numerical values, with the exception of C 4 , 
which in isnr seems to be missing the numerical factor contained in the normalization of He{Z).) The 
coefficients ci. 2 , 3,5 are coupling constants due to volume fluctuations induced by the lowest glueball that 
were apparently dropped in [50]. Coefficients with a star indicate the corresponding constants for the first 
excited exotic mode. 


A. Glueball decay to two pions 


The effective, 3-|-l-dimensional interaction Lagrangian for the lowest (exotic) 0"’“'*' glueball Ge 
reads (omitting terms that vanish when Ge is on-shell) 


_ _rjY 


1 1 

-cid^irduTT^^GE + -cid^ndf^nGE 

E 


(4.2) 


with coupling constants ci and ci defined in (B3) and numerically given in Table III 


The corresponding result for the dilatonic scalar 0'’'+ and the 2'^~^ mode, denoted Ge> and 
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vertex 

value 

di 

Gdirdn 

17.226 A-5 

C ?2 

Gpp 

4.3714 A-5 

da 

Gdpdp 

18.873 A-5 

ds 

Gpirdn 

512.20 A-i Ac"^Mkk 

d) 

G'*i97rd7r 

11.906 A-5 

d 2 

G*pp 

-0.9415 A-5 

dl 

G*dpdp 

13.680 A-5 

dt 

G* pirdiT 

419.46 A-i 


TABLE IV. Coupling coefficients di {ti = V^di) in the interaction Lagrangian of the lowest glueballs in 
the tensor multiplet (dilaton and tensor), collectively denoted as G, with a star indicating the first excited 
mode. (Note that there is no term analogous to the one involving C 4 for the lowest (exotic) glueball.) 


respectively, is 


^ d.nd^TT ( - 


~Mi 

= ltiTrdf,Trd^TTT>^\ ti = V6di. 

Gd is a canonically normalized real scalar, and a massive tensor field with transverse traceless 
polarizations, normalized such that 


Gd, 

(4.3) 

di. 

(4.4) 


+ ... 


(4.5) 


where Bn and B are Lagrange multiplier fields. The coefficient di is given in Table IV 


For the two scalar glueballs described by Gg and Gd, the decay width into two pions is given 
by the simple expression 


Tg 




Jg^|A<E.Dl=x3xl. 


(4,6) 


where p is the momentum of one of the pions in the rest frame of the glueball with |p| = Me,d/‘^, 
the factor of 3 comes from the sum over the isospin quantum number, and the factor of | is included 
because the two pions are identical. The amplitude for the decay of G^; and Go is, respectively. 


^ Mp 

\Me\ = |(ci + ci)poqo - cip • q| = |ci + 2ci|^, 


\Md\ = |diP • q| = |di| 


Ml 


D 


(4.7) 

(4.8) 


For the tensor glueball an average over the polarizations of the tensor is needed. Alternatively, 
we can choose a fixed polarization = 1 and integrate over the orientation of our Cartesian 

coordinates. This leads to the scattering amplitude (in the rest frame of the tensor glueball) 


\-^t\ = \ti{pl - pI)\, |p|= Mr/2, 


(4.9) 


and the decay width 


TT—^TTTT 


SttMI 


47 r' ' 2 6 407r' ' ^ 


(4.10) 
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M 

V/M 

Ge 27r 

855 

0.092 .. 

0.122 

G*^ 27r 

2168 

0.149 .. 

0.197 

Gd 2tt 

1487 

0.009 .. 

0.012 

G*jj 2 tt 

2358 

0.011 .. 

0.014 

T ^2tt 

1487 

0.0145 .. 

0.0193 

T* 2-k 

2358 

0.0175 .. 

0.0233 


TABLE V. Decay width of scalar and tensor glueballs into 2 (massless) pions divided by glueball mass for 
A = 16.63... 12.55 


Numerically we obtairj^with A « 16.63 for the scalar glueballs Ge, Gd, and G*j^ 




3|ci + 2ci|2m| 13.79 


5127r 


AiVi 


Tg 


£)^7r7r 


/Md = 


^G*^^1T7T 


/Md* = 


and for the tensor 






JMt* = 


3|di|2M2 ^ 
5127r 

1.359 
AiV2 ^ 

1.633 

s 0.009 

Ri 0.011 

5127r 

AiV2 


2.174 

! 0.0145 

6407r 

AVI 



2.613 

« 0.0175 

6407r 

AVI 



0.092 {Me ~ 855MeV) 


{Me « 1487MeV) 


{Md* « 2358MeV) 


{Mt ~ 1487MeV) 


{Mt* ~ 2358MeV) 


(4.11) 


(4.12) 


(4.13) 


(4.14) 


(4.15) 


If we replace the standard choice A ~ 16.63 by the smaller value 12.55 as discussed above, all 
these decay rates which are proportional to A“^ increase by 33% (see Table [v| for a summary). 

A somewhat anomalous feature of the lowest (exotic) scalar glueball is that its width is much 
larger than the next-to-lowest (dilatonic) scalar glueball while having a rather low mass. This 
appears rather unnatural if the dilatonic scalar glueball is interpreted as an excited scalar glueball 
and may be another indication that the exotic mode should be discarded altogether. 

Interestingly enough, a scenario with a broad glueball around 1 GeV in combination with a 
narrow glueball in the range predicted by quenched (as well as unquenched 0) lattice gauge theory 
has been proposed in Ref. [25l [75H77] on the basis of QCD spectral sum rules. There the lighter 
glueball, called cjb, plays the role of an important bare glueball component of the cr-meson /o(500), 
while a higher narrow glueball around 1.5-1.6 GeV is required by the consistency of subtracted 
and unsubtracted sum rules. The glueball state a-Q of Ref. [25l [75l [77] has a broad decay width 


into two pions, in fact even much broader than (4.11), which makes us speculate that the exotic 


scalar glueball of the Witten-Sakai-Sugimoto model could hnd a role as the holographic dual of 
a pure-glue component of the cr-meson, perhaps while having to be discarded from the spectrum 
of the pure pure-glue Witten model7*^ This would also be in line with the fact that the gluon 
condensate of the Witten model, Eq. (3.10), is small, close to its standard SVZ value [53j, while 


Ignoring the contribution involving ci, the result for the relative width of the scalar glueball G would read 0.040 
in agreement with the result of [50|, because the fact that the coefficient in |cip is twice that of |50| is exactly 
compensated by A“^ in |ci|^ being half that in [50]. 

This dichotomy might be due to the fact that the flavor D8 branes of the Sakai-Sugimoto model are localized in the 
x'^ direction along which the graviton mode associated with Ge is polarized, whereas this extra spatial direction 
should play no active role in the Witten model—while the requirement of even a;‘*-parity does not rule out the 
exotic mode Ge, some further projection may be appropriate for the pure-glue case. 
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models with only one scalar glueball field |24l l52] cannot reconcile a small gluon condensate with 
narrow glueball states. 

In the range 1.5-1.8 GeV, where lattice gauge theory locates the lowest scalar glueball, there are, 
experimentally, two isoscalar mesons /o(1500) and /o(1710) which are frequently and alternatingly 
considered as predominantly glue. The experimental results for the decay width into two pions are 

r(®^)(/o(1500) ^ 7r7r)/(1505MeV) = 0.025(3), (4.16) 

r(-)(/o(1710) ^ 7r7r)/(1722MeV) = | 

where the first result is taken from Ref. |14) . the second from Ref. m using data from the BES 
collaboration [79] (upper entry) and the WA102 collaboration [80] (lower entry), respectively. 

The lowest (exotic) scalar glueball mode Ge appears to have a much too large decay width to be 
consistent with a dominantly glueball interpretation of either /o(1500) or /o(1710). On the other 
hand, the dilatonic mode has a decay width below but comparable to the data for the two glueball 
candidates; in the case of the WA102 data for the /o(1710) there happens to be even complete 
agreement. In order to get a more complete picture, we shall now consider also the other couplings 
between glueballs and mesons as determined by the Witten-Sakai-Sugimoto model. 


B. Glueball decay to four and more pions 


To leading order in 1/a' or equivalently inverse ’t Hooft coupling, the D8 brane action (3.2) 
does not give direct couplings of glueballs to more than two pions. These appear only through 
higher DBI corrections with terms quartic in the field strength as will be discussed further 
below. 

Decays into more than four pions can however proceed through vertices involving vector mesons. 
The vertices coupling a single glueball to vr and/or p mesons that are obtained from the Yang-Mills 


part of the D8 brane action (3.2) arise from terms of the form (dropping derivatives and Lorentz 
indices) 


GTr(7r7r), GTi{pp), GTr (p[7r, vr]), GTr ([vr, p]^), GTr(p[p,p]), GTr([p,p]^). (4.18) 


Only the first three couplings are relevant for the decay of a glueball to < 4 pions. The correspond¬ 
ing interaction Lagrangians for the exotic and the dilatonic scalar glueball are given explicitly in 
Appendix [B] with the coupling constants for the lowest glueball states listed in Table III and IV 


The relative width of the decay of a glueball to two pions was found above to be Tg^-w-w/M oc 
parametrically suppressed by a factor 1/Nc compared to the decay of the p meson. For 
glueballs with mass larger than 2 mp, the decay into two p mesons is of the same parametric order. 
However, both the lowest exotic glueball and the lowest dilatonic glueball have mass below the 
2p threshold. In this case at least one p meson has to be off-shell, which leads to an additional 
suppression by a factor Tp/nip oc 

Because the vertex coupling a single p meson to two pions involves Tr (p[7r, vrj), the leading- 
order decay into four pions produces pairs of pions with different isospin index. (The parametrically 
suppressed decay G —)• 2G —>■ and G —)• G-|-27r*^ —>■ which only needs the leading Yang-Mills 

part of the DBI action will be discussed together with the direct decay G —?■ Att^ from higher-order 
DBI corrections further below.) 
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FIG. 3. Leading-order glueball decay into four pions, isospin indices a ^b. 


1 . 


Leading-order deeay rate of scalar glueballs to four pions involving ir 


± 


The Feynman diagrams for the amplitude of the decay Ge,d 27r“ -|- with a ^ b are shown 
in Fig. 1^ Some details of the rather lengthy calculation of the decay rate are given in Appendix 

o 


Because one internal p meson can reach its mass shell, while it has nonnegligible width, we 
include (following Ref. [50]) Fp in the p meson propagator according to ^p{r) = l/(rg — — 

rUp -|- impTp) with Fp given by (3.20). This corresponds to a partial summation of higher-order 
terms in inverse powers of XNc- As a crosscheck of our calculations, we have verified that in the 
limit XNc — )• oo the resulting decay rate agrees with the rate for Ge,d pnir, and in the case of 
glueballs above the 2p threshold, with G ^ 2p (Appendix C lb). 


Because rrip < mE,D < ‘2mp, the leading parametric order of the decay width of and Gpi 
into four pions is given by the process G —>• pTrvr and reads X~‘^N~^. Decays through off-shell p 
mesons contribute terms of order X~‘^N~‘^. 


For Ge, which is only 10% heavier than a p meson, the contribution from one on-shell p meson 
is strongly suppressed by phase space, but the finite width of the p meson helps to increase the 
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rate. For A ps 16.63 we fincF^ 

TcE^An/ME ~ 1.33 X 10“^ (A « 16.63). (4.19) 

For the heavier dilatonic glueball Go, the process G —)• pinr is more dominant, leading to a 
significantly larger relative width 

Tco^in/MD ^2A4:XW-^ (A « 16.63). (4.20) 


Evidently, the 47r decay of the lowest holographic glueball state, be it or Gd, is strongly sup¬ 
pressed. Table VI summarizes these results and also shows them for smaller A = 12.55. (In Section 
|V]we shall consider the extrapolation of these lowest states to the higher masses of experimental 
glueball candidates in the range predicted by lattice gauge theory.) 


2. Decay of excited scalar glueballs to two vector mesons 


For the excited dilatonic glueball with mass Mu* ~ 2358.4 MeV, which is above the 2p threshold. 


a similar calculation, but with coefficients dt in place of di (see Table IV), gives 




0.104 (A « 16.63). 


(4.21) 


This result, which involves resummed p propagators, is in fact well approximated by the decay rate 
to two on-shell p mesons: 


14 330 

rG^^^2plMD* « ^ 0.096 (A 16.63), 


(4.22) 


which corresponds to the strictly leading-order part of (4.21) as explained in Appendix C lb 


The result (4.22), divided by its isospin factor of 3, also gives the decay into two isosinglet 


vector mesons w, whose mass is only 1% higher than that of the p meson. 


Since the decay width into two pions given in (4.13) is much smaller than the width into two 


vector mesons, the excited dilatonic glueball turns out to decay predominantly into four pions and 
six pions. 

The excited exotic scalar (if we do not discard this mode altogether) is instead dominated by 
the decay into two pions, which makes this state extremely broad. Calculating also the decay into 
two vector mesons, we find that the decay into two p mesons accounts for only about a third of 
the total decay into four pions. 


^G%^2p/Me* 


2.078 


0.014 (A « 16.63). 


(4.23) 


This means that the decay into four pions is coming largely from the G*^p'K'k vertex. 


In Table VII the results for the decay widths of the excited exotic and dilatonic scalar glueballs 
in the Witten-Sakai-Sugimoto model are summarized. While the excited dilatonic scalar glueball 
has a more moderate decay width compared to the very broad excited exotic scalar, it turns out 
to be still quite large, around 500 MeV. 


Omitting the contributions from the interaction terms involving the c coefficients as in |50j would give the even 
lower value 5.1 x 10“®. In contrast to the decay into 2 pions, in the 4-pion decay rates the factors \/2 in the 
coupling constant qym and in the normalization of the lowest scalar glueball by which we differ from Ref. |50| no 
longer cancel. However, even when using exactly the couplings of Ref. [50] we have not been able to reproduce the 
numerical result 2.2 x 10“® given in Eq. (3.26) of | 5U| . 
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M 

T/M 

> 47r 

855 

1.3 X 10-4 .. 

. 3.0 X 10-4 

Gd — >■ 47r 

1487 

2.4 X 10-3 .. 

. 3.9 X 10-3 

Gd 471° (NLO-DBI) 

1487 

4.0 X 10 -® .. 

. 2.9 X 10-3 

Gd^Ge + 27r° 477° 

1487 

2.6 X 10 -® .. 

. 4.5 X 10-6 

Gd-^Gd+ 277° ^ 477° 

1487 

1.9 X 10 -® .. 

. 4.5 X 10-9 


TABLE VI. Decay widths of lowest exotic and lowest dilatonic scalar glueballs into four (massless) pions 
divided by glueball mass for A = 16.63 ... 12.55. 



M 

T/M 

G%^{27r,2K,2r,} 

2168 

0.397.. 

.0.526 

G*e 47r 

2168 

0.037.. 

.0.061 

G’fj —?► 2u; —>■ Ott 

2168 

0.005.. 

.0.006 

G*e 2(j> 

2168 

0.005.. 

.0.006 

G% (total) 

2168 

0.443.. 

.0.599 

G*jj 47r 

2358 

0.104.. 

.0.142 

— y 2 cj — y Gtt 

2358 

0.032.. 

.0.043 

G*jj 2(1) 

2358 

0.032.. 

. 0.043 

G*j, ^ {2TT,2K,2r,} 

2358 

0.029.. 

.0.039 

GJ) (total) 

2358 

0.197.. 

.0.267 


TABLE VII. Decay widths of excited scalar glueballs divided by glueball mass for A = 16.63... 12.55 (chiral 
limit, with a ratio 3:4:1 for the combined decay into 27r, 2K, 2ri). 


3. Scalar glueball decay to four tt^ 


The glueball decays into four pions that we have considered above involve pairs of pions with 
different isospin index. A decay to four tt^ is suppressed by powers of inverse ’t Hooft coupling, 
because it either has to come from higher-order contributions in the DBI action of the D8 branes 
(Fig. 1^ or has to involve glueball self-interactions and virtual glueballs (Fig. [^. 

As shown in Appendix]^ the parametric order of the vertex formed by a single glueball and four 
TT^ turns out to be whereas the amplitude for G —>• 2 G —)■ and G ^ G + 27r° —)■ 

is proportional to . The former thus has stronger suppression in inverse powers of A, 

while the latter is more strongly suppressed with respect to inverse powers of Nc- 

For simplicity, we only consider the dilatonic glueball, since the exotic glueball has a much more 


complicated interaction Lagrangian. In Appendix B 2 b the interaction Lagrangian for a dilatonic 
glueball with four resulting from the next-to-leading terms of the DBI action has been obtained, 
and in Appendix 


B2c 


the vertex for Gd ^ Gd,e + 27r°. Numerically evaluating the respective 
decay rates of the dilatonic glueball shows that at finite’t Hooft coupling and Nc = 3 the dominant 
decay process comes from the direct coupling of to four vr*^. For A 16.63 we find (see Appendix 


C 2 for details) 


..(NLO-DBI) 

Go^>47r0 


/Md ^ 4.02 X 10“® (A pa 16.63). 


(4.24) 


The decay through virtual glueballs, while not as strongly suppressed by inverse powers of A, 
is subleading at large Nc and is disfavored by phase space. To check whether it might nevertheless 
be important at A^c = 3 and our range of’t Hooft coupling, we have evaluated the first diagram 
in Fig. [^involving one virtual glueball and found that its contribution is smaller than (4.24) by 
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■s. ^ ^ 


Tl 


Tl 


TT 


TT 


FIG. 4. Glueball decay into four 7r° through a vertex from the next-to-leading order terms of the DBI action; 



TT 


TT 



TT 


TT 


TT 


TT 


FIG. 5. Glueball decay into four 7r° (a) through terms in the Yang-Mills part of the DBI action that are 
quadratic in the glueball mode; (b) through a pair of virtual glueballs. 


several orders of magnitude (see Table [V 

^(LO-DBI) 


GD->-GD-l-27rO->-47rO 


/Md ^ 1.94 X 10“^ (A « 16.63). 


(4.25) 


If we do not discard the exotic glueball as a physical state (for instance if we were to interpret 
the latter as holographic dual of a glueball component of the cj-meson, as speculated at the end 
of Section IVA), we should also consider the process Gd —t Ge + 2'7r° which is less suppressed 
kinematically (but still by N~^). This would be of similar magnitude as the result (4.24):^ 

.(LO-DBI) ~ 2.56 X 10 "® 


Gu-S-GB+27r0^.47r0 


/Md 


(A PS 16.63). 


(4.26) 


(As shown in Table VI, at smaller A this contribution is less important compared to the next-to- 


leading DBI contribution (4.24).) 

Decays into four have been seen for the glueball candidate /o(1500) at a level of about an 
order of magnitude below the general dvr decay |14j . whereas no such data seem to be available 


for /o(1710). The smallness of the holographic result (4.24) however would correspond to a much 
stronger suppression than the one observed experimentally for /o(1500). 


4. Tensor glueball decay to two vector mesons 

Unless the mass of the lowest tensor glueball is manually adjusted (as we shall consider to do 
in Section 0 . only the excited tensor glueball of the Witten model with mass Mt* = Md* ^ 
2358.4 MeV can decay into two p or two ui mesons. 

The decay rate involves two sums over the polarizations of the two vector mesons. The average 
over the polarization of the tensor can again be performed by choosing the particular polarization 


By contrast, in the scenario of Ref. [25] . where the cr-meson has a large glue contribution, the heavier glueball is 
claimed to have important decays. 
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= 1 and averaging over spatial directions. The rate for two vector mesons with fixed 
isospin quantum number reads 


r = 


1 


IGttMt’. 


3 

Y{M3,/2)2-m2l'M.x(V,A2)|t 


Ai,A2 = 1 


(4.27) 


where Ai ^2 are labels for the polarizations of the two vector mesons and refers to the specific 
tensor polarization. The amplitude Afe.r(Ai,A 2 ) and the final result of the summations and the 


the result for the decay into two p mesons is 


integration are given in Appendix C 3 With coupling constants t* = and d* from Table 


IV 


Tt *—>.2^—>'47r 

Adj'* 


21.236 

~jNf 


0.142 (A PS 16.63), 


(4.28) 


and 1/3 of that result for the decay T* —)• 2u; —)• Gvr. This should be compared to the decay rate 
into two pions, Eq. (4.15), which is less than 1/8 of (4.28). 

As we shall discuss below, a similar pattern arises when the lowest tensor glueball is extrapolated 
in mass such that it is above the 2p threshold. 


V. EXTRAPOLATIONS AND COMPARISON WITH EXPERIMENTAL DATA 


When comparing our results for decay rates with experiment, it seems reasonable to do so with 
the dimensionless ratio T/M when extrapolating the mass M of the holographic glueball to the 
mass of the experimental glueball candidates /o(1500) or /o(1710). In the case of decay into two 


massless pions, Eqs. (4.11 )-(4.15), this ratio involves two explicit powers of the glueball mass M 


that cancel the inverse mass scale squared coming from the normalization of the glueball field, 
Eq. (2.17) or (2.23). When extrapolating to higher glueball masses, we thus assume that the 


normalization of the glueball field scales according to the glueball mass. While this keeps T/M for 
two-pion decays unchanged, the decay rates into two vector mesons or four pions are modified and 
depend in fact strongly on whether the glueball mass is above or below the 2p threshold. 


A. Extrapolations for the scalar glueball candidates /o(1500) and /o(1710) 


The results of such an extrapolation to the experimental masses of the isoscalar mesons /o(1500) 


or /o(1710) is given in Table VIIT where the holographic results of the (chiral) Witten-Sakai- 
Sugimoto model for the lowest (“exotic”) and the dilatonic 0'*““*' glueball are compared to the 
experimental results for the total and the partial decay widths. Here we have generalized our 
results to Nf = 3 and assumed that pions, kaons and p mesons appear in ratios 3:4:1, respecting 
SU(3) flavor symmetry. 

Explicit masses for quarks would require a modification of the Sakai-Sugimoto model, for exam¬ 
ple along the lines of Ref. [72H71!, which we intend to study in future work. This will necessarily 
modify the coupling of scalar glueballs through contributions that depend on the mass of the 
pseudo-Goldstone bosons, and this may either increase or decrease the decay amplitudes into the 
heavier pseudo-Goldstone bosons. A significant enhancement would be in line with the so-called 
chiral suppression of scalar glueball decays that is suggested by the lattice results of Ref. [HI] and 
the analysis of Ref. [82| . (In the dilaton effective theory of Ref. [52] also an increase of the ampli¬ 
tude for the decay into a pair of heavier pseudo-Goldstone bosons was found, however such that it 
is approximately canceled by the kinematical suppression from the phase space integral.) 
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decay 

^exp 

r/M (exp.) 

T/M[Ge 

(M®^P)] 

T/M[Gd 


/o(1500) (total) 

1505 

0.072(5) 

0.249. 

..0.332 

0.027. 

..0.037 

/o(1500) ^ in 

1505 

0.036(3) 

0.003. 

.. 0.006 

0.003. 

..0.005 

/o(1500) ^ 2^ 

1505 

0.025(2) 

0.092. 

..0.122 

0.009. 

..0.012 

/o(1500) ^ 2K 

1505 

0.006(1) 

0.123. 

..0.163 

0.012. 

..0.016 

/o(1500) ^ 2rj 

1505 

0.004(1) 

0.031. 

..0.041 

0.003. 

..0.004 

/o(1710) (total) 

1722 

0.078(4) 

0.252. 

..0.336 

0.059. 

..0.076 

/o(1710) ^ 2K 

1722 

* f 0.041(20) 

\ 0.047(17) 

0.123. 

..0.163 

0.012. 

..0.016 

/o(1710) ^ 277 

1722 

* f 0.020(10) 

\ 0.022(11) 

0.031. 

..0.041 

0.003. 

..0.004 

/o(1710) ^ 2^ 

1722 

>1= f 0.017(4) 

\ 0.009(2) 

0.092. 

..0.122 

0.009. 

..0.012 

/o(1710) ^ in 

1722 

? 

0.006. 

.. 0.010 

0.021 

..0.030 

/o(1710) ^ 2a; ^ 67 r 

1722 

seen 

0.00016... 

0.00021 

0.011. 

.. 0.014 


TABLE VIII. Experimental data for the decay rates of the isoscalar mesons /o(I500) and /o(17I0) juxta¬ 
posed to the holographic results for the various decay channels of the lowest (exotic) glueball {Ge) and 
predominantly dilatonic glueball (Go) with mass mE,D artificially raised to the respective experimental 
values (still in the chiral limit, i.e. with massless pions, kaons, and ry) and ’t Hooft coupling varied from 
16.63 to 12.55. Experimental data are from Ref. [Hj except for those marked by a star, which are from 
Ref. [75] where the total width of /o(1710) was split under the assumption of a negligible branching ratio 
to four pions, using data from BES |75| (upper entry) and WA102 |5D| (lower entry), respectively. (Holo¬ 
graphic predictions that are substantially increased due to the manually adjusted glueball mass are rendered 
in italics.) 


When comparing the extrapolated decay rates of the holographic glueballs with those of the 
isoscalar mesons /o(1500) or /o(1710) we find that the lowest (exotic) glueball is much too broad to 
be identified as their dominant glueball component. The dilatonic glueball, however, is sufficiently 
narrow for this purpose. It leads to a total decay width that is quite close to the experimental 
width of /o(1710), while being somewhat more strongly below that of /o(1500). With mass equal 
to that of /o(1500), the dilatonic glueball has significantly smaller width in 27r decays, and still 
smaller for decays into four pions, which is the dominant decay mode of the /o(1500). 

Regarding the /o(1710), the decay into 27r is found to be nicely comparable to the experimental 
value, while the stronger rate into pairs of heavier pseudo-Goldstone bosons remains unaccounted 
for with our assumption of SU(3) invariance. A significant enhancement of decays into kaons and 
r] mesons may however be brought about by mass terms for the latter which inevitably will give 
additional contributions to the coupling with scalar glueballs. 

If our extrapolation of the decay width into dvr can be trusted, this appears now uncomfort¬ 
ably large considering that the decay of /o(17IO) into dvr has not been observed. It should be 
noted, however, that the experimental data for the branching ratios of the /o(1710) still have large 
uncertainties and are not covered by the Particle Data Group [Hj. The quoted results are from 
Refs. 1211 [78], which assume that decays into tttt, rjrj, and KK add up to the total width with 
negligible contribution from dvr decays. 

Our extrapolations also predict decays into two cu mesons at a nonnegligible level. According 
to [H], decays of /o(1710) to two u mesons have at least been seen. The Witten-Sakai-Sugimoto 
model for a (pure) glueball candidate suggests that the rate into four pions should be about twice 
as large. 

In this context it is worth mentioning that there are still many open questions surrounding the 
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decay 

M 

r/M[r(M)] 

T^2tt 

1487 

0.013. 

.0.018 

T^2K 

1487 

0.004. 

.0.006 

T ^ 2r] 

1487 

0.0005.. 

0.0007 

T (total) 

1487 

« 0.02 

. .0.03 

T —>■ 2/9 —> 47r 

2000 

0.135. 

.0.178 

T 2K* 2{Ktt) 

2000 

0.119. 

.0.177 

T —>■ 2uj 6 tt 

2000 

0.045. 

.0.059 

T^2tt 

2000 

0.014. 

.0.018 

T^2K 

2000 

0.010. 

.0.013 

T ^ 2rj 

2000 

0.0018.. 

0.0024 

T (total) 

2000 

« 0.32 

. .0.45 

T 2K* 2{Ktt) 

2400 

0.173. 

.0.250 

T —>■ 2p —> 47r 

2400 

0.159. 

.0.211 

T ^ 2uj ^ Qtt 

2400 

0.053. 

.0.070 

T^2(j) 

2400 

0.032. 

.0.051 

T^2tt 

2400 

0.014. 

.0.019 

T^2K 

2400 

0.012. 

.0.016 

T ^ 2p 

2400 

0.0025.. 

0.0034 

T (total) 

2400 

« 0.45 

. .0.62 


TABLE IX. Extrapolation of tensor glueball decay for the case of massive pseudo-Goldstone bosons with 
glueball mass M = Mx = and when the latter is raised to 2 GeV or the lattice prediction ^ 2.4 GeV. 
The’t Hooft coupling is again varied from 16.63 to 12.55. 


nature of /o(1710) [9]. For example, some authors have argued that the nearby resonance |83] 
/o(1790), which is not yet covered by the Particle Data Group, should be combined with /o(1710) 
into one object /o(1760), for which Ref. [83j was able to fit disparate decay patterns, with and 
without significant decay into four pions. 


B. Extrapolations for the tensor glueball 


In the Witten model, with Mkk = 949 MeV, the mass of the tensor glueball equals the mass 
of the dilatonic scalar glueball, and the tensor glueball has roughly similar decay rates into two 
and four pions. The rate into two pions practically exhausts the decays into pions, and has been 


calculated above in Eq. (4.14). The lowest tensor glueball thus turns out to be a rather narrow 


state, however this is due to the fact that it stays below the 2p threshold. 

Indeed, the situation is markedly different for the excited tensor glueball T*. Its mass equals 
that of the excited dilatonic glueball, and because this is above the threshold for two p mesons, 
there is a significant contribution to four-pion decays, and also from other vector meson decays, 
as we have seen in Section IV B 4 Extrapolating the couplings of the lowest tensor glueball to a 


similarly high mass, 2 or 2.4 GeV (where the latter is roughly the prediction of lattice gauge theory 
for the lowest tensor mode), equally gives large contributions from decay into two vector mesons, 
as listed in Table |IX[ Reassuringly, these results are quite close to those for the unmodified results 
for T*, cp. Eq. (4.28), so that we consider them as plausible extrapolations to the likely situation 


of a tensor glueball with mass above 2 GeV. 


In Table IX we have also extrapolated to decays into kaons and p mesons. In the holographic 
setup, a tensor glueball presumably does not couple to an explicit mass term of the pseudoscalar 
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mesons, so the effect of the latter should be purely kinematic. The results (4.9) and (4.10) imply 
a pseudoscalar mass dependence of the form (1 — 4m^/M|.)^/^. This suppression is such that it 
overcompensates the ratio 4/3 that favors kaons over pions. 

For the decays into vector mesons K* and cj) we have taken into account that their masses 
are larger than rup in the phase space factor, but we have left open the possibility that this also 
increases the coupling t 2 = \/ 6 d 2 and merged the two alternatives in the range of results for the 
corresponding decay rates for tensor glueballs with increased mass. 

Adding up the individual contributions, we find a very broad width for a 2.4 GeV tensor glueball, 
1.1 to 1.5 GeV, which is much broader than all the /2 mesons listed in mi. With a mass around 2 
GeV, the width (600 to 900 MeV) turns out to be larger but perhaps marginally comparable with 
that of the tensor meson /2(1950), which has T = 472(18) MeV. The latter is indeed occasionally 
discussed as a candidate for a tensor glueball as it appears to have largely flavor-blind decay modes. 


VI. CONCLUSION 


Using the Witten-Sakai-Sugimoto model for holographic QGD, which only has one free di¬ 
mensionless parameter, we have repeated and extended the calculation of glueball decay rates of 
Ref. [50], where only the lowest scalar mode was studied. 

This lowest mode is associated with an exotic polarization of the gravitational field, involving 
components in the direction of compactification from a 5-dimensional super-Yang-Mills theory 
down to nonsupersymmetric Yang-Mills theory. The mass of this lowest (“exotic”) O"*"^ mode 
turns out to be only slightly above the mass of the p meson and is therefore much smaller than 
the mass scale of glueballs found in lattice gauge theory. The background of the Witten model 
also contains another tower of scalar glueball modes which are predominantly dilatonic and whose 
lowest mass is about 1.5 GeV, not far from the predictions of lattice simulations. 

Besides its very low mass, the lowest (exotic) scalar glueball turns out to have a decay rate 
that is significantly higher than that of the heavier dilatonic mode, which seems counterintuitive 
if the latter were to represent an excitation of the former. We are therefore led to the conjecture 
that the exotic scalar mode should be discarded so that the glueball spectrum begins with the 
(predominantly) dilatonic mode as lowest glueball. Another, more speculative possibility that 
we have mentioned in Section IV A is that the exotic scalar mode represents a broad glueball 

[TSIETI, which features a broad 


component of the cr-meson in line with the scenario of Ref. 
glueball around 1 GeV and a narrower one around 1.5 GeV. 

The decay widths of glueballs obtained in the Witten-Sakai-Sugimoto model are parametrically 
suppressed by a factor of but the numerical results vary substantially for the different 

modes and decay channels, and thus do not give a picture of “universal narrowness” despite the 
large-Vc nature of the Witten-Sakai-Sugimoto model. 

A very strong parametric suppression is obtained for the decay into dvr®, as already pointed out 
in Ref. [5(1] . We have confirmed that also the final numerical value turns out to be very small. 

A noteworthy feature of the Witten-Sakai-Sugimoto model is that the value of the gluon con¬ 
densate is small, close to its standard SVZ value [53] . whereas phenomenological models which 
incorporate a scalar glueball through a QGD dilaton field |24L l52] would require very large gluon 
condensates to admit only narrow glueball states. 

We have also extrapolated our results so that they can be compared with experimental data 
for the scalar glueball candidates /o(1500) or /o(1710). In the case of /o(1500), our results for 
the decay widths of the dilatonic glueball are significantly below the observed rates for decay into 
two pions and even more so for the experimentally dominant decay into four pions. In the case 
of the /o(1710) meson, the decay rate into two pions comes out in nice agreement with available 
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experimental data. The much stronger rate into kaons is not accounted for, but this may be due 
to the fact that the Witten-Sakai-Sugimoto model is strictly chiral and the mechanism of chiral 
suppression [sniHS]. However our (crude) extrapolation to the mass of /o(1710) predicts also a 
significant branching ratio into four pions that has not been seen experimentally. [Although in this 
context it should be noted that the identification of /o(1710) and its separation from the nearby 
/o(1790) [83] has been a matter of debate P l8l]]. 

Furthermore, we have studied the decay of tensor glueballs, which in the Witten-Sakai-Sugimoto 
model have a narrow width into two pions and (when the mass is above the 2p threshold) a large 
width into four pseudoscalars, such that at best the isoscalar tensor meson /2(1950) appears to be 
(marginally) compatible with our the holographic result, while heavier tensor glueballs would have 
to be broader than the tensor mesons so far discussed in the literature. 

In the case of the tensor glueball we can already plausibly anticipate the effects of nonzero 
pseudo-Goldstone masses. In the case of scalar glueballs the situation is less clear and we intend 
to study this issue in extensions of the Witten-Sakai-Sugimoto model in a future work. This would 
be particularly interesting in view of the glueball candidate /o(1710) which according to Ref. |21| 
could be a nearly unmixed glueball and which has a ratio [Hj r( 7 r 7 r)/r(iFiF) that is significantly 
below the flavor-symmetric value 3/4. 

Since the holographic results pertain only to pure glueballs, it would clearly be most interesting 
to study mixing of glueballs with qq states as this can strongly obscure signatures of glueball 
content. In the holographic setup, mixing is suppressed by 1/W |SQ| and would presumably 
require more difficult stringy corrections that are not captured by the effective Lagrangian following 
from the Witten-Sakai-Sugimoto model. Absent those, it might be interesting to consider a more 
phenomenological approach such as extended linear sigma models [ 23 ], where holographic results 
for the glueball-meson interactions could be used as input instead of fitting to experimental data. 
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Appendix A: Ten-dimensional field equations 


The Kaluza-Klein reduction of the eleven-dimensional graviton modes yields metric fluctuations 
pertaining to the four-sphere, i.e. in the components denoted by Omitting these fluctuations 
as done in Ref. [50] corresponds to dropping all vertices proportional to c* in the interaction 


Lagrangian (B2) of the exotic scalar glueball. To see whether this reduction could be justified, we 


check if these truncated modes solve the ten-dimensional field equations of type IIA supergravity. 

For the Witten-Sakai-Sugimoto model, the relevant terms of the supergravity action are given 
by |85| 


5iia = ^ j (^R + 4Vm , (Al) 

where 2K^g = (27r)^Zf and F 4 = dC^ is the four-form from the R-R sector of the theory, with 


iFif = j^FabcdF^^^^. 


(A2) 
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Variation of this action with respect to the background metric qmn and the dilaton field $ 
results in 


1 e^® 

Rmn - -^Rqmn + 2VM<hVAr<l>- — FmABC “ 3!|iV|^) = 0 (A3) 


and 


R + 4VmV^$ - 4Vm = 0, 


(A4) 


respectively. 

The solution to these equations that corresponds to the background of the Witten model is 
given by the metric (2.7), the dilaton (2.8) and a nonvanishing R-R four-form held. The latter is 
hxed by the requirement that the hux through a unit four-sphere is quantized, i.e. 


[ F^ = N,^, (A5) 

JS4 va TT 

where the factor of Nf. arises from the fact that we are considering a stack of W D4-branes. The 
held strength that satishes this condition is given by F 4 = 3i?D4<?7^a;4, with a ;4 denoting the volume 
form of the unit four-sphere 

With this information, one can linearize the held equations and plug in the solutions both with 
and without the spherical huctuations, which is easily done with computer algebra tools. The 
result is that for both the dilatonic and exotic glueball modes, the held equations are not satished 
unless the huctuations along the four-sphere are included. This means that in a rigorous top-down 
approach the vertices corresponding to the coefficients c* have to be included in the calculation 
of decay rates. (For the dilatonic glueball mode, the need to include huctuations can also be 
deduced from the explicit 10-dimensional calculations of Ref. [39j.i 


Appendix B: Glueball-meson interaction Lagrangians 


The effective interaction Lagrangian of glueballs and qq mesons is obtained by inserting the 
10-dimensional metric huctuations (4.1) into the D8 brane action and integrating over the bulk 
coordinates. In this section we give the result for the lowest (exotic) scalar glueball, the dilatonic 
scalar glueball, and the tensor glueball, expanded up to the order needed for the calculation of 
decay rates of glueballs into two pions and four pions as discussed in the text. 

As discussed above, we do so by keeping induced huctuations in g^Q. In the D8 brane action 
(3.1) the contribution from the dilaton huctuation dGnqi appearing through the factor ^/gs^ = gfiQ 
in y/—g is opposite in sign to that from e~^ and larger by a factor 4/3. 

Let us also recall that following Ref. |50j we use the convention 


7r = 7r“r“, /9^ = p“T“, Tr^“^^ = (5“^ (Bl) 

so that for Nf = 2 using Pauli matrices we have r“ = and Tt\T°-,T^]T^ = while 

Ref. |48( I49| have TrT^T^ = The Minkowski metric used in the 3-|-1-dimensional Lagrangians 
is = diag(- -k -k-k). 


Note that this result looks different in some of the literature, e.g. |48]. This is due to a different convention with 
rescaled three-form potential. 
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1. Lowest scalar mode 

The glueball-meson interactions contributing at leading order to the decay of a glueball into two 
or four pions are given by the terms linear in the glueball field and up to quadratic in p, maximally 
trilinear in tt, p in the Yang-Mills part of the DBI action of the D8 branes. In the case of the lowest 
(exotic) scalar mode, they read 


^Ge = -Tr < Cl 


df^nd^TT^^GE + \{df,7rf Ge 

E ^E-' 


+C2 M^k 


■'e 




+C 3 


Ip 

^ ^ UO^ V 






+icr> 


d^d'' 

9.vr[vr,p.]^ 


Ge + ^9;,7r[7r,p^] (^1 - Ge 


+ ^ci5^7r(9'"7r Ge + ^C2 mIj^p^p^^Ge 
+ \c2,F^yF^''GE + ic59/,7r[7r,p^]GE >, 


(B2) 


where F^y = dfj_py — d^p^j, without a commutator term This agrees with Ref. [50], whose 

notations we have adopted, in the part involving c,, but Ref. [50] effectively dropped all terms 
proportional to Ci due to the neglect of dgo,^. 

The coefficients Cj and q are obtained by integrals over the glueball mode function He{Z), the 
p meson mode function ppZ), and the pion mode function (j)Q{Z) oc 1/K, K = 1 + Z'^, according 
to 


Cl = / dZ^, C2 = 2 kJ dZK{p[fHE, C3 = 2 kJ dZ K-^^^{PifHE, 

O T\j 2 f 20ZK , f PiHe 

C 4 = 2 kMkk J dZj^^^^—^pip^HE, J ’ 

Cl = / dZ^^, j dZK{p[fHE, 

C3 = l^^j dZK-^^\PifHE, h = j dz"^^, (B3) 


where the integral over Z is from 


—00 to -|-oo and where following Ref. |50| we have introduced 


He{Z) = 


1 , 3 

4 5R:- 2 


He{Z). 


(B4) 


The corresponding coefficients for the excited mode G*^ are obtained by replacing the lowest 
mode function He{Z) by the next highest eigenfunction. 

The numerical results for the coefficients Ci, Ci for the lowest mode as well as for c\ and for 
GT are given in Table III 
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2. Dilatonic and tensor mode 


a. Gluehall-meson interactions contributing to leading order decays 


Restricting ourselves again to glueball -meson interactions contributing to leading order decays 
to two and four pions, the interaction Lagrangian linear in Gd or T, up to quadratic in p, and 
maximally trilinear in vr, p reads, for the dilatonic mode, 


j Gd + d2M-^^-p^j,py f - — 


Cgd = Tr di-df.'Kdy'K ( - 


Gd 


1 - - / d^d^\ ( 

+d^-F,pF,P - — j GD + id,d^7r[7r,p,] - — ] Gd 


(B5) 


with coefficients 


d^ = jdZ^, d2 = 2K j dZK{'ilj[fHD, 
d3 = 2Kj dZK-^/\i^ifHD, d5 = I dz'f^, 

and for the tensor glueball 

£t = Tr I d^nd^TT 

+ \h F^pF.PT^^'' + ih 9 ^ 7 r[ 7 r, p,] 


(B6) 


(B7) 


with ti defined in analogy to (B6). Because of Ht oc. Hd with the normalization conditions (2.23), 
we simply have ti = \/6di. 

The numerical results for d* and the corresponding coefficients d* for the next-highest dilatonic 


scalar are given in Table IV 


b. Gd-4t^'^ vertex from next-to-leading order DBI action 

A direct coupling of glueball modes to more than two pions appears only at higher orders of the 
DBI action of the D8 branes. For the coupling to four vr^ = vr^ we need to expand up to quartic 
terms in Fj^^- The action, restricted to Fj^^, reads 

5 = rD8(2W)2 I • (B8) 

Inserting the metric fluctuations corresponding to the dilatonic glueball and dropping terms 
that vanish on the mass shell of the glueball gives 

^ 3^/ [{d^f^GD - 2d[idyfidy)id,7r^) Gd (B9) 

with 

4 = 8A3^Jm4^ / dZRDK-^!^ ^ 2.513 • 10® A'i (BIO) 
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c. Two-glueball-two-TT vertices 


The leading (Yang-Mills) part of the DBI action also contains nonlinear terms with respect to 
the metric fluctuations dual to glueballs, which have to be considered for the glueball decays in 
four TT^, which vanish at leading order. Expanding the bilinear term in ir to second order in the 
dilatonic mode Gd yields 

■■GoGoTrTr _ r 




-TV 


3{d^7r)^Gl 


-d^irdyir rip^ - 


□ 


dudf 

Goir''-^] Gd 


with 


df^ = 


dZ^ = 399.04... A"^ 
ttK 


(Bll) 


(B12) 


In Eq. (4.26) we have considered for completeness also the decay through the lowest exotic 
scalar glueball. Eor this process the relevant terms in the interaction Lagrangian turn out to be 


^GdGbttOttO ^ c^^Tr 


dpird’^TT 


1 




--di.T^d^TT 


( Vi 

fdPd^^ 


Ge I ( Vo-u 




□ 


Gd 


Vcrp 




(B13) 


with 


cf ® = 1653.9 ... A'^iVT^Mr^ 


ttK 


KK- 


(B14) 


Appendix C: Four-pion decay amplitudes and phase space integrals 
1. Decay of scalar glueballs into 4 massless pions involving 


The leading-order decay amplitude of a glueball into four pions involves two pairs of pions with 
different isospin index (thus excluding the case of four vr’^’s). If Ai is the amplitude for G —i 27r“27r^ 
with fixed a ^ b, the total decay rate of a glueball into 4 pions is given by 

rG^4. = 1x^1 dLlPS,{M)\M\^ (Cl) 

where the factor | is due to a factor of 3 for the three different pairs a, b possible, and | is the 
symmetry factor for two pairs of identical particles. 

Eor the decay of a particle at rest with mass M into n particles we have 

rfLIPS„(A/) = (mss - E !>a) n 


Useful details of how to organize the integration over the final momenta are given in Ref. [50] . 
As a test of the numerical procedure for implementing these integrations we have used that for 
massless final states the phase space integral with A4 = 1 can be done analytically with the result 

m 


j dLIPSn(M) 


j\^2n-4 

2(47r)2n-3r(7j)r(n — 1) 


(C3) 
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a. Four-pion decay amplitude for the dilatonic glueball 

Because the lowest glueball corresponding to an “exotic” polarization of the metric fluctuations 
has a rather lengthy interaction Lagrangian and because we arrived at the conjecture that the next- 
higher scalar (predominantly dilatonic) mode should be interpreted as the lowest scalar glueball of 
QCD, we shall give the decay amplitude into 4 pions explicitly only for the latter. Denoting the 
final pion four-momenta in Gd -f 27r“27r^ by p,p', q, q' (see Fig. and defining 

= = r^=p'^ + g'^, s^^=p^^ + qf^, r° + (C4) 

we find 


i.AA. — ^s/^QpTVK ) “1“ ^5 ^ * t) 

+9%n^pir)Ap{s)\ d2Ml^a • b 


+^3 


{aobo — a-b)r-s-|-a-b (roso — r • s) 

— (aoro — a • r) b • s — (feo'So — b • s) a • r 

+ {q q') ) 


(C5) 


where Ap(r) = l/(rQ — -|- irripYp) with Tp given by (3.20). 


b. Scalar glueball decay through pnir and pp 

The use of a finite width of the p meson in the propagator Ap corresponds to a partial resumma¬ 
tion of formally higher order diagrams. This seems to be natural in view of the fact that Tp/rrip is 
not very small, but it should be kept in mind that e.g. a correction of the residue of the propagator 
is being dropped. 

If the glueball decay were to be treated strictly perturbatively in inverse powers of A and 
one would neglect Ypjmp as a higher-order contribution and treat the p meson as nearly stable. 
Because to leading order there is no local vertex that would couple a glueball directly to four 
pions, the leading-order process would then be given by a decay into on-shell pinr with decay 
width proportional to as long as the glueball mass is below the 2 p threshold; glueballs 

with mass larger than 2 mp would have the decay into two p mesons as dominant process for the 
eventual decay into four pions, with partial width proportional to . 

We have evaluated the decay rates into pinr (and pp when M > 2mp) as a cross-check of our 
results for the decay into four pions, which coincide in the limit of large A, 

lim X^Yg^Att/M = lim X^Yg^p-k-k/M = 71 for rrip < M < 2mp, (C 6 ) 

A—^00 A—^00 

lim ArG'-s. 47 r/Tf = lim XYg^ 2 p/M = 72 for M > 2mp, (C7) 

A—^00 A—^00 

with 7 i ~ 22.074iV“^ for the glueball mode Gp), and 72 ~ 6.451Ai“^ when its mass is artificially 
raised to 1722 MeV. Taking these strictly leading-order results as a basis for the decay width into 
four pions would give somewhat higher numerical values than the above calculation involving a 
finite Tpj^ 


Ref. |50] has added FG-^ 47 ^ and Fwhen comparing their results with experimental data, which we regard as 
overcounting. 
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For the lowest (exotic) glueball mode, whose mass is not much higher than rup, we find 71 ~ 
0.0030iV“^ (~ 0.00114A^“^ if the contribution from the c’s is dropped). Again the results for 
converge to this limit for A —)• 00 . However, for A = 16.63 the effect of resumming 
Vp in the calculation of rG'_^ 4 ^ is now an increase of the decay width by more than two orders 
of magnitude compared to the strictly perturbative result that corresponds to a nearly stable p 
meson with negligible width: the latter would give rG^^p-Kn^iTr/ME ~ 4.0 x 10“^ compared to 
1.3 X 10“^ from Tge^Ati/Me with resummed p propagators. 

The decay amplitudes for G —)• pTrvr are somewhat unwieldy, in particular for the exotic glueball 
mode. We therefore give details only for the decay of the excited dilatonic glueball into two p 
mesons. (For the analogous decay of the lowest dilatonic glueball when its mass is raised above 
the 2 p threshold is obtained by replacing d* by di.) 

No phase space integration is involved in this process, but the polarizations of the p meson have 
to be summed over. Denoting the two transverse and the one longitudinal polarization by indices 
T and L, respectively, the result is 

Te*^ppIMd* = (|M^|2 + 2\Mt\^) , (C 8 ) 


with 


I-^lI = 


d*.M? 


KK 


Ml. 

4mp 


+ dim. 


\Mi 


^2-^KK + *^3 



(C9) 


Decays into two u) mesons, whose mass equals the p meson mass in the Sakai-Sugimoto model 
(in the real world it is only 1 % heavier), is given by the same expression with the overall isospin 
multiplicity factor 3 omitted. (For the excited dilatonic glueball, also decay into two 4> mesons 
becomes relevant, though not for the lowest dilatonic glueball when its mass is raised to one of the 
glueball candidates as the latter are all below the 2 (f) threshold.) 


2. Decay of dilatonic glueball into 4 massless 7r° 


With M. the amplitude for the total decay rate is given by 


rG,^4.o = ^ X ^ /dLIPS4(M)|M|2. 


The dominant contribution is provided by (B9), which leads to 


M. / {Sid'i) = 3 [p ■ qp ■ q + p ■ p q ■ q + p ■ q' q ■ p] 
—p • qp' • qf — p ■ qp' ■ q' — p ■ p' q- q' 
-p ■ p' q ■ q' - p ■ q' q - p' - p ■ q ■ p'■ 


(CIO) 


(Cll) 


3. Decay of tensor glueball into two vector mesons 


Unless one adjusts its mass parameter, only the excited tensor mode is above the 2p threshold. 


The interaction Lagrangian (B7) with U replaced by U determines the amplitude in Eq. (4.27) 


for a specihc tensor polarization ejf and two p mesons with momenta p, q and polarizations 
e^(p, Ai),e''(g, A 2 ) as 


A4ej,(Ai,A2) — ep{p, Xife^iq, X2) 


-tl {p - (T ■ qp^'" +P- qf-T - P'" ^'t <lp - <1^^7 Pp) 


(C12) 
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With = M/2, p = —q, and = (M/2)^ — the result of the summation over the 

polarizations Ai, A 2 and the integration over spatial directions reads 



+4^2^3-^kk (^3P^ 4- ' m - p ^ 

+^tf (;^(P^)^ + 2m^ P^ + ”^p) > 


with t* = V^d* and d* given in Table 


IV 
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